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ABSTRACT 

| We have resimulated the six galaxy-sized haloes of the Aquarius Project including 

metal-dependent cooling, star formation and supernova feedback. This allows us to 
study not only how dark matter haloes respond to galaxy formation, but also how 
this response is affected by details of halo assembly history. In agreement with pre- 
y's !■ vious work, we find baryon condensation to lead to increased dark matter concen- 
tration. Dark matter density profiles differ substantially in shape from halo to halo 
when baryons are included, but in all cases the velocity dispersion decreases monoton- 
ically with radius. Some haloes show an approximately constant dark matter velocity 
| anisotropy with f3 » 0.1 — 02, while others retain the anisotropy structure of their 

baryon-free versions. Most of our haloes become approximately oblate in their inner 
regions, although a few retain the shape of their dissipationless counterparts. Pseudo- 
phase-space densities are described by a power law in radius of altered slope when 
baryons are included. The shape and concentration of the dark matter density profiles 
are not well reproduced by published adiabatic contraction models. The significant 
, spread we find in the density and kinematic structure of our haloes appears related to 

■ differences in their formation histories. Such differences already affect the final struc- 

ture in baryon-free simulations, but they are reinforced by the inclusion of baryons, 
and new features are produced. The details of galaxy formation need to be better un- 
derstood before the inner dark matter structure of galaxies can be used to constrain 
| cosmological models or the nature of dark matter. 
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pie, the slow increase of rotation velocity with radius at the 

centre of low surface brightness galaxies (Flores & Primack Galaxy formation might reinforce such history-specific 

1994; McGaugh & de Blok 1998; Moore et al 1999; Salucci, features in the dark matter distribution, since baryons are 
Yegorova & Drory 2008) and the relatively weak concentra- subject to dissipative processes such as cooling, star forma- 
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of baryons within dark matter haloes modifies both their 
dynamics and their structure, but exactly how this happens 
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the dark halo is compressed radially and adiabatically by 
the added mass at its centre (Blumenthal et al. 1986, here- 
after B86; see Eggen et al. 1962, Zeldovich et al. 1980 and 
Barnes & White 1984 for previous applications of this for- 
malism) . Recent simulations have found the very simple adi- 
abatic compression (AC) scheme of B86 to overestimate the 
effect (e.g. Gnedin et al. 2004 and Abadi et al. 2009, here- 
after G04 and A09, respectively; Pedrosa, Tissera & Scan- 
napieco 2010) but it is nevertheless often used in simplified 
modelling of galaxy formation (e.g. Mo, Mao & White 1998). 
More sophisticated AC models, based, for example, on the 
formalism of Young (1980) have claimed to reproduce better 
the contraction and the final shape of haloes (G04; Sellwood 
& McGaugh 2005). The common outcome of these schemes 
is nevertheless an increase in the dark matter density in 
the central regions, exacerbating the problem of reconciling 
ACDM models with observation. 

Simulations of galaxy formation which follow both dark 
matter and baryons in their proper cosmological context 
are the primary tool for studying how galaxy assembly af- 
fects dark matter haloes. There are many papers devoted to 
this subject. Improving numerical algorithms and increas- 
ing computer power have led to continual progress in under- 
standing the complex interplay between these two compo- 
nents. Previous analyses of the evolution of dark haloes in 
hydrodynamical simulations have reported an increase both 
of the central mass concentration and of the central velocity 
dispersion, which no longer shows the 'temperature inver- 
sion' characteristic of dissipationless CDM haloes (e.g. Katz 
& Gunn 1991; Evrard, Summers & Davis 1994; Navarro & 
White 1994; Tissera & Domingucz-Tenreiro 1998; Onorbe 
et al. 2008; Romano-Diaz et al. 2008). Recently, Pedrosa, 
Tissera & Scannapieco (2009, 2010) investigated how bary- 
onic assembly history affects the dark matter distribution. 
They found the SN feedback process to play a key role by 
regulating star formation activity and ejecting material not 
only from the main system, but also from infalling satellites. 
It seems clear that the final distribution of baryons at the 
centre of a halo is insufficient to determine halo response to 
galaxy assembly, thus contradicting the AC hypothesis. 

Typical dark matter haloes have long been reported to 
be triaxial with major to minor axis ratios often exceed- 
ing two in CDM-only simulations (e.g. Barnes & Efstathiou 
1987; Frenk et al. 1988; Dubinski & Carlberg 1991; Jing 
& Suto 2002; Hayashi, Navarro & Springel 2007). However, 
their shape changes, becoming more nearly oblate as baryons 
condense within them (Katz & Gunn 1991; Evrard et al. 
1994; Tissera & Dommguez-Teneiro 1998; Kazantzidis et al 
2004; Debattista et al. 2008; A09). These results suggest the 
need for a comprehensive analysis of halo structure in sim- 
ulations which include both the appropriate cosmological 
context and a realistic description of the physics of baryon 
condensation. It is clearly necessary to analyse a number of 
different galaxy-sized haloes in order to explore how differing 
assembly histories affect the structure of the final haloes. 

In this paper, we study a set of high-resolution resimula- 
tions of the six galaxy-mass haloes of the Aquarius Project 
(Springel et al. 2008). These were carried out with a ver- 
sion of GADGET-3 which includes a multi-phase treatment 
of metal-dependent cooling, star formation and SN feedback 
(Scannapieco et al. 2005, 2006). The original haloes were se- 
lected from a cosmological CDM-only simulation with no re- 



striction on merger history, except that implied by eliminat- 
ing objects with high-mass close neighbours. This Aquarius 
halo set is well-suited to study how formation history and 
baryonic condensation together determine the final structure 
of dark matter haloes. Our ability to isolate these effects is 
aided by comparing results from our hydrodynamical simu- 
lations (hereafter SPH runs) with corresponding results from 
the CDM-only simulations (hereafter DM runs) as reported 
by N08. Properties of the galaxies in these haloes are stud- 
ied in Scannapieco et al. (2009) and will be further analysed 
in a forthcoming paper. 

This paper is organized as follows. In Section 2, we de- 
scribe the numerical experiments. In Section 3, we analyse 
the dark matter density profiles. Section 4, describes the ve- 
locity dispersion structure of our haloes. In Section 5, we 
study the effects of baryons on the pseudo-phase-space den- 
sity profile, while Section 6 discusses halo shapes. Section 
7 compares the change in halo circular velocity profile be- 
tween the SPH and the DM runs with the predictions of 
AC models. Finally, in Section 8, we summarize our main 
findings. 



2 THE SIMULATED HALOES 

The six haloes studied in this paper were taken from the 
Aquarius Project (Springel et al. 2008). They were se- 
lected at random from a lower resolution version of the 
(100 /i _1 Mpc) 3 Millennium-II Simulation (Boylan-Kolchin 
et al. 2009) subject only to the requirements that their mass 
should be similar to that inferred for the Milky Way's halo 
and that they should have no close massive neighbour. Dark 
matter only versions of these haloes were then resimulated 
at a variety of much higher resolutions as part of the Aquar- 
ius Project itself (Springel et al. 2008; N08; Vogelsberger et 
al. 2009). The lowest resolution version of each halo (des- 
ignated as resolution level 5) was also resimulated with de- 
tailed modelling of baryonic processes by Scannapieco et al. 
(2009, hereafter S09). Details of how the initial conditions 
were created can be found in these two papers. In the follow- 
ing we will refer to the original Aquarius simulations at the 
higher resolution level 2 as the DM simulations, and to the 
simulations of S09, which have roughly 200 times worse dark 
matter mass resolution but include the baryonic physics, as 
the SPH simulations. The analysis in N08 shows that the 
properties of the dark haloes which concern us in this paper 
are extremely well numerically converged all the way down 
to resolution level 5. The individual haloes are labelled Aq- 
A-n, Aq-B-n, etc., with n corresponding to the resolution 
level in to be consistent with the convention in the earlier 
papers. 

The Millennium and Millennium-II Simulations, and 
thus also these halo simulations, were carried out assum- 
ing a ACDM cosmology with parameters f2 m = 0.25, J1a = 
0.75,o-8 = 0.9, n s = 1 and H = 100 h kms^Mpc -1 with 
h = 0.73. Evolution in both the DM and the SPH simula- 
tions was followed from z = 127 to z — using versions of 
GADGET-3, an update of GADGET-2 (Springel et al. 2001; 
Springel 2005) optimized for massively parallel simulation 
of highly inhomogeneous systems such as individual dark 
haloes. 

The version of GADGET-3 used for the SPH runs 
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includes a multiphase model for the gas component 
with metal-dependent cooling, star formation and phase- 
dependent treatments of SN feedback and chemical enrich- 
ment, as set out in Scannapieco et al. (2005, 2006) and em- 
ployed for studying ACDM galaxy formation in Scannapieco 
et al. (2008, 2009). The model describes chemical enrichment 
by SNII and SNIa separately, using appropriate yields and 
delays (Mosconi et al. 2001; Scannapieco et al. 2005). The 
multiphase treatment of both hydrodynamics and feedback 
is quite effective in reproducing the observed phenomenol- 
ogy of star formation and wind-generation in both quiescent 
and starburst galaxies. Without introducing system-specific 
parameters, the scheme produces substantial mass-loaded 
galactic winds in rapidly star-forming systems with speeds 
that reflect the escape velocity much as in observed galaxies 
(Scannapieco et al. 2006, 2008). This freedom from ad hoc 
scale-dependent parameters, makes the algorithm particu- 
larly suited for studying galaxy formation in its cosmologi- 
cal context, since the simultaneous formation of systems of 
widely differing mass is the norm in this situation. 

The SPH runs have maximum gravitational softenings 
in the range cg = 0.5 — 1 ft _1 kpc. To prevent spurious re- 
sults due to limited numerical resolution in the innermost 
regions, we analyse halo properties only outside 2/i -1 kpc 
in all our systems. Haloes are considered bounded by their 
virial radius r2oo defined as the largest radius within which 
the mean enclosed density exceeds ~ 200 times the critical 
density. The simulated SPH haloes have s=s 1 million parti- 
cles in total within this virial radius while the DM simula- 
tions have approximately 10 s particles in the same region. 
The virial masses of the six systems are in the range 5 to 
11 x 10 11 /i _1 M Q . Hence, dark matter particles in the SPH 
runs have masses of the order 10 6 /i _1 Mq while gas particles 
initially have « 2x 1O 5 /i _1 M . In Table [Qwe summarize the 
principal characteristics of our SPH simulations. We also in- 
clude the corresponding information for two lower resolution 
versions of halo Aq-E. We find convergent results between 
Aq-E-5 and Aq-E-6 for their characteristic properties. How- 
ever, for the lowest resolution run, Aq-E-7, larger differences 
are found. In Section 3, we discuss the effects of numerical 
resolution on the dark matter distributions. 

The Aquarius haloes have varied assembly histories and 
these produce a variety of structures and star formation his- 
tories for their central galaxies, even though the simulations 
were all run with similar SF and SN feedback parameters. A 
detailed description of the code and the parameters adopted 
can be found in S09, together with images of the central 
galaxies and considerable analysis of their structure. Here 
we provide a brief summary so that the reader has a general 
picture of the central galaxies of the haloes we analyse here. 

In all six galaxies star formation peaked between 13 Gyr 
(Aq-C-5) and 10 Gyr (Aq-B-5) ago and there has been rather 
little star formation over the last 8 Gyr. The total stellar 
masses at z = are in the range 2.5 to 6 x lO lo /i _1 M0 and 
are listed in Table [1] along with the total masses and baryon 
masses of the z = haloes and the numbers of simulation 
particles representing these masses. 

The analysis of S09 showed that most of the cen- 
tral galaxies contain centrifugally supported disks, although 
none of these accounts for more than about a fifth of the 
total stellar mass. Aq-F-5 is an exception in that it has no 
disk, only a spheroidal component. This can be traced to the 



fact that the system underwent a major merger at z ~ 0.6. 
In some cases the stellar spheroid has substantial nett rota- 
tion (Aq-E-5) while in others it retains very little angular 
momentum (Aq-A-5). All the galaxies have been substan- 
tially affected by SN feedback, which drives winds which 
limit their mass. In all cases the baryon fraction within the 
virial radius is around 10%, significantly smaller than the 
global value, /b = 0.16. 



3 DARK MATTER DENSITY PROFILES 

We have constructed spherically-averaged dark matter den- 
sity profiles from 2 /z -1 kpc to the virial radius r2oo- The in- 
ner radius is four times the gravitational softening for three 
of our simulations (A, B and D) and twice the gravitational 
softening in the other three. It typically contains more than 
a thousand dark matter particles. Finding an appropriate 
centre is crucial to obtaining accurate profiles. Here we use 
the shrinking sphere technique as laid out in of Power et al. 
(2003). In practise, this gives a centre very close to the min- 
imum of the gravitational potential in all our systems. We 
have measured profiles both including and excluding sub- 
structures, but this only causes rather minor effects in the 
outer regions. In the following we use the full DM mass 
including substructures in order to be able to compare con- 
sistently with the results of N08. Note that we re-normalize 
the dark matter profiles of N08 by the global DM fraction 
of 0.84 adopted in the SPH simulations when comparing 
profiles from the two set-ups. 

To characterise our dark matter profiles we adopt the 
three-parameter Einasto model (Einasto 1965) which N08 
found to give a relatively good fit to the profiles of the DM 
runs. This model has also been used by earlier authors to 
fit the dark matter distribution when baryons are included 
(e.g. Onorbe et al. 2008; Gao et al. 2008). N08 parametrise 
the Einasto profile using a, r_2 and p_2, which indicate its 
curvature in a log-log plot, and the radius and density at the 
point where its logarithmic slope is —2, the isothermal value. 
We fit this formula, leaving all three parameters free, to our 
measured density profile at a set of points spaced logarith- 
mically from 2 ft _1 kpc to r2oo, minimising the rms residual 
in log p. In Table[2]we give the results of these fits, together 
with the rms residual in log p. To facilitate the comparison 
with the DM runs, we also include the corresponding param- 
eters given by N08; note that these were obtained for fits over 
a slightly wider radial range. We see that in all cases r— 2 
and a are both smaller for the simulations which include 
baryons, showing that the condensation of the galaxy has 
increased the concentration of the halo and given it a den- 
sity profile which is approximately isothermal over a wider 
radial range than in the DM-only case. The residuals in log p 
are only slightly larger, showing that the Einasto profile is 
a good fit to both the SPH and the DM runs. 

This is evident in Fig. [1] where we compare the 
spherically-averaged dark matter density profiles of the 
haloes in our SPH runs (solid lines) to the corresponding 
profiles for the DM runs from N08 (dashed lines). Every 
SPH profile lies below its corresponding DM profile at large 
radii and above it at small radii. In all cases the two profiles 
cross at, or slightly inside, the baryonic radius of the central 
galaxy, which we define to be the radius containing 83% of 
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Table 1. General characteristics of our Aquarius halo simulations and their central galaxies. The first column show the halo name, taken 
from Springel et al. (2008). ec ls the gravitational softening. r2oo and M200 are the radius and mass of the sphere enclosing total mean 
mass density 200 times the critical value. Afdm and Mt, ar are dark matter and baryonic masses within this same sphere. M s is the stellar 
mass of the main galaxy hosted by each halo, measured within an optical radius defined to contain 83% of the galaxy's baryonic mass. 
iVdm and iVbar are the total number of particle representing each mass component. 
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Table 2. Characteristics of the density profiles of the haloes in our Aquarius galaxy formation simulations. Column (1) gives the name 
of each halo. Columns (2), (3) and (4) list a, r_2 and p_2 (in units of kpc h~ 1 and MQh 2 kpc — 3 ), the parameters of the best fitting 
Einasto model. Column (5) gives the rms scatter in logp around this fit. Column (6) lists exponents for the best power law fits to p/<r 3 
as a function of r. Columns (7), (8) (9) (10) and (11) list the corresponding fitting parameters for the corresponding DM runs. Note that 
the fitting parameters for the profiles in the DM runs differ slightly from those given by N08 because we have re-calculated them over 
the same radial range used for the SPH runs. 
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its stars and cold gas. Inside the baryonic radius, the SPH 
profiles are all relatively flat, and so approximately "isother- 
mal" . This is indicated more explicitly in the insets in Fig. [1] 
which plot the logarithmic derivatives of the two profiles as a 
function of radius. This slope changes smoothly with radius 
in most of the DM haloes, but shows a clear change in be- 
haviour in the inner regions in many of the SPH haloes. For 
Aq-A-5, Aq-C-5 and Aq-E-5, we detect a change in the log- 
arithmic derivatives of the SPH profiles at the smallest radii 
suggesting a flattening of the profiles. However, this occurs 
very close to our resolution limit and therefore, should be 
viewed with caution. At radii well outside the central galaxy 
the SPH profiles parallel the DM profiles quite closely, and 
indeed the offset between them is close to the renormalisa- 
tion factor that we introduced to account for the difference 
in the total amount of dark matter in the two kinds of sim- 
ulation. 

We have looked for correlations between the changes in 
a and r-2 and various other parameters of our haloes, in 
particular, their spin parameter, the fraction of the baryons 
gathered in the central galaxy, and their disk-to-total mass 
ratio. However, in no case did we find a clear trend. 

To study the robustness of our results against numerical 
resolution, we analysed the dark matter profiles of the two 
lower resolution versions of halo Aq-E. The results show ex- 



cellent convergence for the dark matter profile and its char- 
acteristic properties, as can be seen in Fig. \T\ Runs Aq-E-5 
and Aq-E-6 are in very good agreement over the whole anal- 
ysed region. The same is true for the lowest resolution Aq-E- 
7 but only over the region resolved with at least a thousand 
dark matter particles. A lower resolution seems to produce 
an underestimation of the dark matter density in the inner- 
most regions. In Table 2, we include the fitting parameters 
for these two additional simulations (carried out from twice 
their gravitational softening). Based on these findings, we 
are confident that our main set of simulations have been 
run with adequate numerical resolution to assure conver- 
gent results for the density profiles over the analysed radial 
interval. 



4 VELOCITY DISPERSION DISTRIBUTIONS 

The velocity dispersion structure of the dark matter haloes 
in the DM runs has been analysed by N08, who show that 
they all have a temperature inversion in the central region 
(i.e. the velocity dispersion drops towards the centre at small 
radii) but that there is no simple regularity among the 
haloes. Rather there is considerable diversity which seems 
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Figure 1. Spherically-averaged dark matter profiles for six haloes from the Aquarius Project simulated with (solid lines) and without 
(long dashed lines) baryons. We plot the profiles from 2 h~ 1 kpc to the virial radius r200- The DM runs are the resolution level 2 
simulations studied by Navarro et al. (2008). The arrow indicates the baryonic radius of the central galaxy defined as the radius enclosing 
83% of its stars and cold gas. The inset boxes show the logarithmic derivatives of the profiles. We include the lower resolution versions 
of Aq-E: Aq-E-6 (short-dashed lines) and Aq-E-7 (dotted lines). 
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to reflect the variety of formation paths possible for these 
galaxy-sized haloes in the A-CDM cosmology. 

We here carry out a similar analysis of the velocity dis- 
persion structure of our set of SPH runs. In Fig. [2] we com- 
pare their total velocity dispersion profiles to those of the 
DM runs. It is clear that the baryons have substantially af- 
fected the velocities of dark matter particles within the bary- 
onic radius of the central galaxy. As reported in previous 
work (e.g. Tissera & Dominguez-Tenreiro 1998; Romano- 
Diaz et al. 2008; Pedrosa et al. 2010) the increased velocity 
dispersion in the central regions results in a final profile that 
decreases monotonically with radius, in contrast to the be- 
haviour seen in DM-only simulations. For Aq-E-5, we have 
included its lower resolution versions (Aq-E-6 and Aq-E-7). 
As can be seen, the lowest resolution version provides a poor 
representation of the velocity distribution obtained with the 
higher resolution runs. However, Aq-E-6 and Aq-E-5 agree 
very well except in the very central region. 

The conventional velocity anisotropy parameter, de- 
fined as fi(r) = 1 — of/(2of) where a t and ay are tan- 
gential and radial dispersions averaged over spherical shells, 
gives some indication of the internal orbital structure of the 
haloes. As shown by N08, the DM runs are almost isotropic 
(P ~ 0) in their inner regions, then in most cases become 
progressively more radially biased (/? > 0) with increasing 
radius out to about r_2 before becoming more isotropic 
again at even larger radii. There are however, substantial 
differences between the individual haloes. We find that bary- 
onic effects modify the dispersion structure within these 
haloes in a complicated way, as can be seen in Fig. [3] While 
some of the SPH haloes have similar anisotropy structure to 
their DM counterparts (Aq-A-5, Aq-C-5 and Aq-D-5), oth- 
ers have become less dominated by radial motions (Aq-F-5, 
Aq-B-5 and Aq-E-5). All our haloes remain slightly radi- 
ally biased, /3 ~ 0.1 — 0.2, in the inner resolved regions. 
The lower resolution versions of Eq-E (Aq-E-6 and Aq-E-7) 
show a large level of noise around the trend determined by 
the highest resolution run (Aq-E-5). In order to get further 
insight into the effects of baryons, we plot of against of. in 
Fig.|4]both for the SPH runs and for their DM counterparts. 
For the DM runs, nearly isotropic behaviour is evident in 
the cool central regions (r < r_2; magenta circles, radius in- 
creasing with dispersion) with an inversion of the relation in 
the outer regions where radial motions dominate (blue cir- 
cles radius increasing with decreasing dispersion). Baryon 
condensation (either by infall or mergers) increases both ra- 
dial and tangential dispersions and changes the shapes of 
orbits in the central regions. (Note that the axis ranges of 
the main plots and of their insets are different). However, 
each halo has its own particular distribution. Haloes Aq-B-5, 
Aq-E-5 and Aq-F-5, which have approximately constant ve- 
locity anisotropy, are also those that have an almost isother- 
mal density profile to radii larger than the baryonic radius 
of their central galaxy, as is evident from the insets of Fig. 

m 

Finally, we look at the relation between the logarithmic 
slope of the density profile (7(f)) and the velocity anisotropy 
parameter (f3(r)). Hansen & Moore (2006) suggested that 
there might be a "universal" relation between these pa- 
rameters in dissipationaless haloes. N08 found that their 
DM haloes follow a well-defined relation only in the cen- 
tral regions where both parameters are monotonic functions 



of r. In Fig. [5] we show the anisotropy-slope relation for 
our SPH runs. As can be seen, there are two different be- 
haviours. Three of the haloes (Aq-A-5, Aq-C-5 and Aq-F-5) 
show similar and well-defined relations over the rough range 
2 fe _1 kpc < r < r_2, which are fairly close to the HM06 re- 
lation. However, for r > r_2, each halo behaves differently. 
In general, they show large variations in f3 (which differ from 
halo to halo) but relatively minor variations in 7 (compare 
with Fig. [TJ. For the other haloes (Aq-B-5, Aq-D-5 and Aq- 
E-5) the logarithmic slope shows only a weak dependence on 
radius in the region dominated by baryons, consistent with 
a nearly isothermal profile, while the anisotropy varies from 
~ 0.1 to ~ 0.4. As a result the HM06 relation fails dramat- 
ically in these cases. 



5 THE PHASE-SPACE DENSITY PROFILE 

Taylor & Navarro (2001) noticed that the quantity Q(r) = 
p(r)/a(r) 3 is well approximated by a power law Q(r) oc r~ a 
with a ~ 1.875 over more than 2.5 orders of magnitude 
in radius in simulated dark matter haloes of any mass in 
a CDM cosmology (see also, Rasia, Tormen & Moscardini 
2004; Ascsibar et al. 2004; Vass et al. 2008, N08). This is 
the behaviour predicted in the spherical similarity solution 
of Bertschinger (1985). This Q(r) is related to the entropy 
distribution within haloes, but the origin of its power law 
behaviour and its relation to structure formation are not yet 
understood. The apparent universality is established only for 
pure dark matter haloes and may be broken by the effects of 
baryons which modify both the velocity dispersion and the 
density profiles of the dark matter. 

As may be seen in Fig. [SJ where we show Q(r) for 
both the SPH and the DM runs, all our SPH haloes show a 
less steep relation than their DM counterparts or predicted 
by Bertschinger's (1985) similarity solution. The residuals 
(small boxes in Fig. [6} show that the modified profile is 
still well fit by a power law at most radii although in a cou- 
ple of cases (in particular Aq-A-5) there is some indication 
for deviations in the innermost region. In Table [2] we give 
the best fit power-law indices for the SPH haloes over the 
radial range 2 /i _1 kpc < r < r^oo- Interestingly, but not 
surprisingly considering our previous findings, the level of 
departure from the Bertschinger value depends on the halo. 

It could be possible that the contraction of the dark 
matter proceeds adiabatically and that the detected increase 
in entropy at each radius merely reflects the contraction it- 
self. To investigate this point, we plot Q as a function of 
enclosed dark matter mass (relative to the total dark mat- 
ter mass within r2oo) rather than as a function of r. If the 
contraction were purely adiabatic, the SPH haloes should 
show identical behaviour to their DM counterparts when 
Q(r) is plotted against M(r)/M2oo in this way. As shown in 
Fig. [7] this is not the case in any of our haloes. 

The lower numerical resolution versions of Aq-E do not 
show significant differences in the pseudo-phase-space den- 
sity profiles with respect to Aq-E-5 as can be seen from Fig. 
[6] because the numerical effects in the density and velocity 
profiles compensate each other. However, if this relation is 
plotted as a function of the enclosed mass as seen in Fig. [7] 
the differences are slightly larger. 
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Figure 2. Velocity dispersion as a function of radius for the SPH (solid lines) and DM (dashed lines) haloes. The arrows indicate the 
baryonic radii of the central galaxies. We include the lower resolution versions of Aq-E: Aq-E-6 (short dashed lines) and Aq-E-7 (dotted 
lines) . 




Figure 3. Anisotropy parameter (/J) as a function of radius for the SPH (solid lines) and DM (dashed lines) runs. The arrows indicate 
the baryonic radii of the central galaxies. We include the lower resolution versions of Aq-E: Aq-E-6 (short dashed lines) and Aq-E-7 
(dotted lines). 
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Aq-A-5 Aq-D-f 




3.8 4.0 4.2 4.4 4.6 4.8 5.0 3.8 4.0 4.2 4.4 4.6 4.8 5.0 

log 2a \ log 2o \ 

Figure 4. Tangential dispersion versus radial dispersion for the SPH haloes (velocities are in units of km s _1 ). The corresponding 
relations for the DM haloes of N08 are shown in magenta for 2 h _1 kpc < r < r_2 and in blue for r_2 < r < r2oo m the insets. The 
solid black line depicts equality. The SPH haloes show a monotonic relation with radius for both components with the higher dispersion 
values for smaller radius as shown in Fig. [2] 
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Aq-D-5 Aq-E-5 Aq-F-5 




Figure 5. Velocity anisotropy /? as a function of the local value 7 of the logarithmic slope of the density profile in the SPH runs. Filled 
symbols correspond the radial range 2 /i _1 kpc < r < r_2 and open circles to the range r_2 < r < T200- The Hansen & Moore (2006) 
relation in shown as a solid line. 





Figure 6. Pseudo-phase-space density profiles as a function of radius for haloes in the SPH (solid lines) and DM (dashed lines) runs ( 
Q is in units of h 2 MQkpc~ 3 (km s~ 1 )~ 3 ). Residuals of the two relations from the power-law relation predicted in Bcrtschingcr's (1985) 
similarity solution are shown in the narrow lower panels of each plot. We include the lower resolution versions of Aq-E: Aq-E-6 (short 
dashed lines) and Aq-E-7 (dotted lines). 
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Figure 7. Pseudo-phase-space density profiles as a function of enclosed dark matter mass fraction for haloes in the SPH (solid lines) 
and DM (dashed lines) runs (Q is in units of h 2 Mokpc — 3 (km s" 1 )" 3 ). The profiles are shown over the same radial range already shown 
in Fig. [6] We include the lower resolution versions of Aq-E: Aq-E-6 (short dashed lines) and Aq-E-7 (dotted lines). 



6 SHAPES 

Haloes formed by hierarchical clustering of dark matter in 
CDM scenarios are typically triaxial (e.g. Barnes & Efs- 
tasthiou 1987; Frenk et al 1988; Jing & Suto 2002) with the 
inner regions being approximately prolate (e.g. Hayashi et 
al. 2007). However, when baryons condense to make galax- 
ies in these inner regions they become rounder and ap- 
proximately oblate (Dubinski 1994; Tissera & Dommguez- 
Tenreiro 1998; Kazantzidis et la. 2004). This latter result 
seems in better agreement with a variety of observational 
shape estimates for galaxy haloes (e.g., Sackett & Sparke 
1990; Kuijken & Tremaine 1994; Koopmans et al. 1998; 
Helmi 2004; Weijmans et al. 2008). We measure the shapes 
of our SPH and DM haloes following the Dubinski & Carl- 
berg (1991) method, based on the eigenvalues of the moment 
of inertia tensor, /. For each bin of N particles we computed 
ellipsoidal radii 



r = I x + 



il/2 



(6/o 



(c/o 



(1) 



The semi-axes of the triaxial ellipsoids (a > b > c) are 
calculated iteratively as 



a In 



in 



(2) 



where In > I22 > I33 are the eigenvalues of the tensor of 
inertia {Ijj, = Y}. x\x\ lr 2 ). In order to determine a, b and 
c, an iterative cycle is set up, starting with b/a = c/a — 1. 

In Fig. [8] we show the axis ratios and the triaxiality 
parameter T defined as T = (a 2 — b 2 )/(a 2 — c 2 ) as a function 



of radius. Purely prolate objects have T = 1, while purely 
oblate ones have T = 0. In agreement with previous work, 
we find a general trend for the SPH runs to be more oblate 
than their DM counterparts. Although in all haloes both axis 
ratios increase at most radii, each halo shows its particular- 
ities. The largest change in shape occurs in Aq-C-5 which 
becomes oblate with T < 0.4 over the whole analysed range. 
The smallest change in shape is found in Aq-F-5 which only 
weakly modifies its overall prolate shape, although the axis 
ratios do increase noticeably. The only galaxy without any 
disc component at all inhabits this halo. Halo Aq-E-5 shows 
a strong dependence of T on radius in the DM run, varying 
from T w 0.6 to T ~ 0.1. This behaviour is preserved when 
baryons are included although the axis ratios do increase 
slightly. Interestingly, the galaxy in this halo is the only 
one with a strongly rotating bulge. Initially triaxial shapes 
become much more nearly oblate in the baryon dominated 
regions for Aq-A-5, Aq-B-5 and Aq-D-5. These haloes have 
galaxies with extended, diffuse surviving discs. For halo Aq- 
E we have included the analysis of the shapes for the two 
lower resolution versions in Fig. [S] In agreement with the 
other parameters, the shapes tend to converge as numerical 
resolution increases. As expected, the dark matter halo in 
the lowest resolution version, Aq-E-7, has the noisiest shape. 



7 CIRCULAR VELOCITY AND THE 

ADIABATIC CONTRACTION MODELS 

The adiabatic contraction (AC) hypothesis has been used 
extensively to model the effects that baryon condensation 
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Figure 8. Axis ratios b/a (solid lines) and c/a (dotted-dashed lines) for haloes in our SPH (black) and DM (green) runs. In the right 
panels, the triaxiality parameter T is also shown; T = for oblate systems while T = 1 for prolate ones. We include the lower resolution 
versions of Aq-E: Aq-E-6 (blue) and Aq-E-7 (red). 



produces on dark matter haloes. As discussed in the Intro- 
duction, there is now a consensus that the simplest imple- 
mentation of this hypothesis overestimates the dark mat- 
ter concentration in the central region and makes it very 
difficult to reconcile observations with CDM models. This 
discrepancy may well reflect issues other than dark matter 
compression by baryonic gravity, but adopting an inappro- 
priate AC hypothesis can certainly exacerbate the problem. 
In this section, we will apply the standard model of B86, the 



modified algorithm of G04 and the statistical motivated re- 
lation of A09, comparing their predictions to our simulated 
dark matter circular velocity profiles. 

The differences between these various recipes can be 
seen in Fig. [9] These plots show the directly measured 
dark matter circular velocities (V c = [GA/d m (r)/r] 1//2 where 
Mdm(r) is the dark matter mass enclosed within r) in our 
SPH (red lines) and DM (black lines) simulations and com- 
pares them to the curves predicted by modifying the result 
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Figure 8 -continued 



for the DM run according to each of the three recipes, as- 
suming the baryon distribution found in the SPH run. The 
DM profiles have been reduced by 8% to account for the 
global baryonic fraction of 0.16 adopted in the SPH simu- 
lations. In the narrow panels at the bottom of each plot we 
show the residuals from the pure DM result for each of the 
models (the three dashed curves) and for the dark matter 
distribution actually found in the SPH simulation (the solid 
red curve). As expected, the AC model of B86 overpredicts 
the amount of dark matter in the central regions by up to 
50% or so. The G04 and A09 models provide a better repre- 
sentation with A09 being closer to the "true" result in most 
cases. However, for two haloes (Aq-B-5 and Aq-D-5), the 



A09 model underestimates the dark matter mass by several 
tens of percent at some radii. This is expected, since the 
Abadi et al (2009) prescription is based on a large sample of 
dark matter haloes, and for individual haloes one can expect 
overestimates and underestimates in roughly equal numbers, 
just as we find. Note that for Aq-B-5 and Aq-D-5 the G04 
prescription actually gives the best description at all but the 
smallest radii. 

In general terms, the dark matter mass in the innermost 
regions is usually overpredicted and, more importantly, the 
shape of the circular velocity is often poorly matched by 
all of these models. At least in three of our haloes (Aq- 
A-5, Aq-C-5 and Aq-E-5), the velocity curves rise in the 
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Figure 9. Dark halo circular velocity curves predicted by applying the AC models of Blumenthal et al. (1986, violet dashed lines), 
Gnedin et al. (2004; blue dashed lines) and Abadi et al. (2009, green dashed lines) to the circular velocity curves of our DM simulations 
assuming the central baryon distribution found in our SPH simulations. We also show the directly measured circular velocity curves in the 
SPH (red solid lines) and DM (black dotted-dashed lines) simulations. Thus if the AC models worked perfectly their curves would match 
the red solid lines. In the lower panels we show residuals of the three model predictions from the original DM curves, and compare them 
with the residuals for the actual curves in the SPH simulations. For this comparison, the curves for the AC models and DM simulations 
have been scaled down by 8% to account for the baryon fraction in the SPH simulations. We also show circular velocity curves for the 
total mass (DM + baryons) in the SPH runs (magenta lines). 
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central region more gently than predicted by the models (in 
agreement with Fig.l). 

8 CONCLUSIONS 

We have analysed the dark matter distributions in six 
galaxy-sized haloes belonging to the Aquarius Project, com- 
paring results from the original dark matter only simulations 
to those from re-simulations including baryonic processes 
(S09). In agreement with previous work, we find that dark 
matter haloes become more concentrated when baryons con- 
dense at their centres, but that the characteristics of the 
contraction do not correlate in a simple way with the total 
amount of baryons. Our main contribution here is to anal- 
yse similar mass haloes with a variety of formation histories 
using high resolution simulations with and without baryons. 
Our results show that the response of haloes to the presence 
of baryons is sensitive to the details of halo assembly. This 
set of six galaxy-sized haloes provides the opportunity to 
study which properties can be considered common to such 
haloes and which depend significantly on their particular 
formation history. 

Our findings can be summarized as follows: 

• In the regions dominated by baryons, haloes become 
significantly more concentrated than their dark matter only 
counterparts. The level of concentration varies significantly 
from object to object, however, it is not simply related to 
the total baryonic mass accumulated in the central galaxy. 
For r < r_2, the dark matter density profiles of many of 
our simulations are nearly isothermal except, possibly, very 
close to the centre. 

• The velocity dispersion structure is modified in all 
haloes, with velocity dispersion increasing monotonically to 
small radii in all cases. The temperature inversion of NFW 
profiles is no longer present once the galaxy has formed. In 
some systems, the tangential dispersion increases more than 
the radial dispersion (but not all), causing them to become 
more nearly isotropic. 

• The /3-7 relation proposed by Hansen & Moore (2006) 
is obeyed at most over a restricted radial range for some 
of our haloes. It works best in the central region of those 
systems where the isothermal behaviour extends over a rel- 
atively small range. The departures from their predictions 
become large when the profile is isothermal or has constant 
anisotropy over an extended radial range. 

• Pseudo-phase-space density no longer follows the same 
power law in radius as in Bertschinger's (1985) similarity 
solution once galaxy formation is included. The profile dif- 
fers from one halo to the next and in no halo is it consistent 
with the purely adiabatic contraction of the dark matter 
only case. 

• As in previous work, the condensation of baryons makes 
the central regions of all our haloes less aspherical and more 
nearly oblate, although in two cases the changes are small. 
One of these has no significant disk component and the other 
has a bulge with substantial nett rotation. 

• None of the simple adiabatic contraction models pro- 
posed in earlier work is able to describe how the radial 
density profiles of our haloes are modified by baryon con- 
densation. The scheme suggested by Abadi et al (2009; see 
also Pedrosa et al. (2010)) is the most successful of those 



we consider, although it can significantly over- or underes- 
timate the effects in individual cases. For the same haloes, 
the central mass is overestimated by more than 50%. The 
circular velocity curves of our galaxy formation simulations 
are not as centrally peaked as predicted by such adiabatic 
contraction models, but are still more rapidly rising than in 
many observed spirals. Since none of our simulated galaxies 
is disk-dominated it is unclear whether this is a problem. 

Our analysis show that haloes of similar virial mass 
respond in different ways to baryon condensation depend- 
ing on the details of their assembly history. Consequently, 
a much more detailed understanding of galaxy formation 
is needed before we can make reliable predictions for the 
detailed structure of the dark matter haloes surrounding 
galaxies. 
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